Twisted-light-induced optical transitions in semiconductors: Free-carrier quantum 

kinetics 
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We theoretically investigate the interband transitions and quantum kinetics induced by light 
carrying orbital angular momentum, or twisted light, in bulk semiconductors. We pose the problem 
in terms of the Heisenberg equations of motion of the electron populations, and inter- and intra- 
band coherences. Our theory extends the free-carrier Semiconductor Bloch Equations to the case of 
photo-excitation by twisted light. The theory is formulated using cylindrical coordinates, which are 
better suited to describe the interaction with twisted light than the usual cartesian coordinates used 
to study regular optical excitation. We solve the equations of motion in the low excitation regime, 
and obtain analytical expressions for the coherences and populations; with these, we calculate the 
orbital angular momentum transferred from the light to the electrons and the paramagnetic and 
diamagnetic electric current densities. 
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I. INTRODUCTION 



in 1992fi 



Since the seminal work by Allen et al. 
there has been a steady increase of interest in the the- 
ory, experiments, and applications of light carrying or- 
bital angular momentum (0AM), or twisted light (TL)j2, 
Studies in this area span several subfields of physics, 
such as research on the classical/quantum properties of 
TLf^"— its generationji and the interaction of TL with 
atoms/molecules^"— and Bose-Einstein condensates 
At the same time, the interaction of general inhomoge- 
neous light beams with solids is becoming an active field 
of research toof^^"— 

Recently, we laid down the basic theoretical elements 
to study the interaction of semiconductors and insula- 
tors with confined beams of twisted light We obtained 
the optical transition matrix elements of the TL-electron 
interaction and studied the transfer of orbital angular 
momentum using a simple, perturbative wave-function 
approach. That approach was adequate as a first theo- 
retical step, but it has a number of limitations. Being a 
single-particle theory, is has the drawback of not taking 
into account the Pauli exclusion in the photo-excitation 
of multiple electrons, and furthermore, it leaves out the 
electron-electron interaction effects. Thus, a more com- 
plete theoretical treatment of the interband excitation by 
twisted light of solids is called for. 

In this paper we develop a set of "twisted-light- 
generalized semiconductor Bloch equations" (TL-SBE) 
from the Heisenberg equations of motion of the popula- 
tions and coherences of the photo-excited electrons. This 
theory is valid for pulsed or CW twisted-light beams, and 
takes fully into account the inhomogeneous profile of the 
beam, as well as the transfer of momentum from the light 
to the electrons in the plane perpendicular to the beam's 
propagation direction. As long as excitonic phenomena 
are not targeted, the Coulomb interaction does not play 



an essential role in the basic physics of band-to-band op- 
tical transitions, and for that reason we will limit our- 
selves, for the time being, to a free-carrier formulation of 
the theory. From a practical point of view, we mention 
that the free-carrier theory is already involved enough to 
merit a separate presentation, obviously 0jS &- first step 
in a program that aims at obtaining and solving, first 
the mean-field TL-SBE, and later the same equations 
with collision termsJ^ Collision terms describe the scat- 
tering processes undergone by the photoexcited electrons, 
namely, electron-electron, electron-phonon, and electron- 
impurity scattering. Collision terms in the relaxation- 
time-approximation can be added straightforwardly to 
our theory in order to describe qualitatively those scat- 
tering processes, and a numerical solution of the resulting 
equations of motion would allow us to explore the infiu- 
ence of collisions on the effects described here. We leave 
this numerical treatment for future work, which will in- 
clude, besides collisions, the study of strong and pulsed 
TL excitation. Finally, notice that while we concentrate 
here on bulk systems, our theory can easily be applied to 
two-dimensional systems excited at normal incidence. 

Usually, the optical excitation in bulk systems is the- 
oretically dealt with by assuming that the system is a 
cube, quantizing the electrons using cartesian coordi- 
nates and taking, at the right moment in the derivation, 
the limit of large system size. For symmetry reasons, this 
method allows straightforward calculations in the case of 
excitation with plane-wave light. However, it leads to 
a cumbersome formulation in the case of excitation by 
twisted light. This is clearly so because the twisted light 
beam has an inherently cylindrical nature. In two previ- 
ous works on the interacion of TL with quantum dots^S 
and quantum ringS)^ the cylindrical nature of the TL 
beams was handled conveniently by also using cylindri- 
cal coordinates in the description of the electronic states. 
In the theory presented here for bulk systems, we take 
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advantage of this simple but key idea. We imagine the 
sohd as a cyhnder, quantize the electron states in cylin- 
drical coordinates, and finally take the limit of large sys- 
tem; we rely on the fact that bulk properties are then 
independent of the geometry of the solid. We keep, nat- 
urally, the required microscopic structure of the Bloch 
wave functions in order to characterize the valence and 
conduction band states: the periodic parts of the Bloch 
states are approximated by their values at zero crystal 
momentum, a common practice known as effective-mass 
approximation. The use of cylindrical rather than carte- 
sian coordinates allowed us to decouple the Heisenberg 
equations of motion according to values of the electron 
angular momentum, which greatly reduces the complex- 
ity of the problem. Using these generalized TL-SBE we 
predict the kinetics of electrons, show the occurrence of 
electric currents with complex profiles, and demonstrate 
the transfer of 0AM from the light beam to the electrons. 

The paper is organized as follows. The TL vector po- 
tential and the system Hamiltonian are given in Sec. HIl 
Section Hill presents the generalized TL-SBE in terms of 
cylindrical electron states, the partial decoupling of the 
equations of motion, and the perturbative solution. The 
electron quantum kinetics is analyzed with the help of 
the electric current and transferred angular momentum 
in Sec. IIVI Conclusions are given in Sec. |Vl 



II. SYSTEM AND TWISTED LIGHT 

In this work we consider a direct-gap semiconductor or 
insulator and study interband transitions caused by illu- 
mination with a beam of twisted light. We assume that 
the light's frequency is such that mainly band-to-band 
transitions occur so that exciton creation is unimportant. 
Under this regime it is satisfactory to formulate our the- 
ory not including the Coulomb interaction between car- 
riers. Thus, our theory can describe the electrons' ki- 
netics, from irradiation to a fraction of picoseconds — to 
avoid strong deviation due to decoherence — in a large 
number of physical systems, from semiconductors having 
band-gap of a fraction of an eV (e.g. InSb with 0.23 eV) 
through several eV (e.g. GaN with 3.5 eV), up to insula- 
tors having larger gaps, provided that the frequency of 
the twisted field is tuned above the energy bandgap; this 
requires the use of twisted fields in the near-infrared to 
UV spectrum, which docs not constitute an experimental 
difficulty. We are particularly interested in characteriz- 
ing the transfer of angular momentum between the TL 
and the electrons, and describing what the electron dis- 
tribution looks like as a result of the photo-excitation. 
Although several valence bands may be involved in inter- 
band optical transitions, here we consider for simplicity 
a two-band model. The generalization of the theory to 
the case with more than one valence band involved is 
straightforward. 

The vector potential of the TL beam in the Coulomb 



jauge is given by2S 

A{r,t) = Aoe^(«^^-"*) [e±Mqrr)e''^T 



ie.^J;±i(<Z.r)e'('±i)'^ 

1z 



+ c.c. , 



(1) 



with the polarization vectors given by e± = x ± iy = 
e^"^(f±i0) and c.c. denoting the complex conjugate. In 
Eq. ([1]), the radial profile of the beam is given, for con- 
creteness, by Bessel functions, Ji{qrr) and Ji±i[qj.r) — 
an alternative formulation would use Laguerre-Gaussian 
modes instead^. 

The light-matter interaction is described using the 
minimal-coupling Hamiltonian, whose dominant contri- 
bution (for moderate field intensities) comes from the 
linear term p • A. If the TL beam is such that qr < Qz, 
which is verified for usual sizes of the beam's waist, the 
largest coupling term comes from the transverse compo- 
nent of A(r, t), and we only need to consider 

A(r,i) ~ e.Aoe^^"^^-"*) J/(g,r)e''* + c.c. 

= AW(r,t) + A(-)(r,t), (2) 

with a = ±. 

The lowest-order light-matter interaction Hamiltonian 

is 

q 



q 



-p. A(r,t) 



AW(r,f) + A(-)(r,t) 

(-) 



(3) 



with p the momentum operator, and rrie and q the mass 
and charge of the electron. (Note that the nie appearing 
in this equation is the bare electron mass and not the 
effective mass.) We emphasize that, unlike most work in 
light-matter interaction, we must keep the spatial varia- 
tion of the field, in order to capture the relevant physics. 
Thus, Eq. encodes all multipolcs, as is clearly seen 
from the Power-Zienau-WooUcy transformation<2^ 

The complete electronic Hamiltonian in second quan- 
tization for the general multi-band case is 



n 



ba 



{b'a'\hj\ba}a 



t 

b' a' 



a-ba 



ba,b' ' 



where b, b' denote energy bands, like heavy hole, light 
hole, conduction, etc., a is a collective index of quantum 
numbers appropriate for the problem at hand, and a^/a 
are creation/annihilation operators. 



III. FREE-CARRIER SEMICONDUCTOR 
BLOCH EQUATIONS 

A. General equations of motion 



al,^,aba- The 



Let us consider the operator pb'a',ba = "■^a'^'ba- -luc 
equation of motion for this operator in the Heisenberg 
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picture is 



ifi-rPb'a'ba = [pb' a' .ba ■ 

at 



(4) 



For concreteness, we take one type of circularly polarized 
light, either a+ or (7_. While, say, a+ light connects both 
light holes and heavy holes to conduction band states, 
these two processes remain independent of each other 
during the evolution under the Hamiltonian that we are 
considering. Then, for circularly polarized light, we can 
accurately describe the electron kinetics under the opti- 
cal excitation within a two band model. We now special- 
ize Eq. (jH) to a two-band case by considering the evolu- 
tion of the three types of operators Pca'.ca, Pva',va, and 
Pva'.ca, where v (c) stands for the chosen valence (con- 
duction) band. After expanding the commutators and 
assuming that the interaction hi connects only valence- 
to conduction-band states, the equations of motion be- 
come 



a' ,ca ^ca.ca' Pea' .ca ~i~ 

Ca\ hi \vai)pca' ,vai - {vai \ hi \ca') Pi,a^^ca, (5) 



at 



— ^voL^va' Pva.' .va ~\~ 

va\ hi \cai) pya' ,cai - {cai\ hi \va')pcai,va, (6) 



i.Pva'ca 

at 



A 



coL^va.' Pva' ,ca 



E' 

ai 



Ca\ hi \vai)pya',vai 



[cai 



hi \va')pcai,ca , (7) 



where Aba,b'a' = £ba - Sb'a'- 

In what follows we work with the equations of motion 
of the expectation values of the operators p: 



Pc,a'a — {Pca' .ca} : 
Pv,a'a — {Pva'.va} ^ 
Pva'.ca — {Pva' ,ca} : 



(8) 



where the average (. . .) is taken over the initial state of 
the material. These expectation values represent pop- 
ulations when they have repeated indices and quantum 
coherences when they are off-diagonal matrix elements. 
Notice that in Eqs. ([5]l7]) we keep the intrahand coher- 
ences, which are essential in the TL excitation process. 
These coherences are usually left out of the theory when 
the vertical transition approximation is made. 



B. Electronic states in cylindrical coordinates 

The semiconductor Bloch equations for bulk systems 
are usually formulated in the basis of electronic Bloch 
states given by 

^6k(r)-(r|6k) = -^e"'--«,k(r), (9) 



where h is the band index, k is the crystal momentum, 
and L is the linear size of the system. In principle, this 
basis set could be used as well in our treatment of TL- 
excited systems, but one finds that it is not a convenient 
choice. To see this, let us recall the interband matrix 
element of the TL-matter interaction Hamiltonian given 
in Ref. Ill: 



(ck'|4+Vk) 



; g Aa{t) 1 

191/ „ \ —i u t 



where k = k' — k has an azimuthal angle 9 and its 
projection in the x-y plane has length k^, and pb'b = 
(l/a^) J^d^xul,{r)pu},{r). After inserting this matrix 
element (and its complex conjugate) into Eqs. ©-([T]), 
one immediately realizes that a single state in the va- 
lence (conduction) band is connected to a multitude of 
states in the conduction (valence) band having any aiigle 
6; pictorially, this has been represented by us in Ref. [l6[ 
by a cone-like excitation in momentum representation. 
Thus, the equations are almost completely coupled — 
with the exception of the z component — and even the 
perturbation-theory solution looks complicated and hard 
to interpret. 

As anticipated in the introduction, there are com- 
pelling reasons to adopt a different basis set for the elec- 
trons in the solid: i) the symmetry exhibited by the vec- 
tor potential Eq. ([3]); ii) (bulk) properties are the same 
for a box- or cylinder-shaped solid, in the large-system 
limit; and Hi) optical excitations are well described in a 
band-edge or effective mass approximation. Therefore, 
we adopt cylindrical states to treat the electrons; their 
wave functions and energies are (see Appendix |X| 



^hk,n(r) = AAj„,(A;,r)e^™'^e'^^^Mfa(r), 

h"^ 

She Eg , 



Sbkr. 



2mt 



(k -\-- k 



(10) 



where k represents the quantum numbers {fc^, fc^}, with 
kz = 2'Kn/H and kr = rms/Li^ H and L are the height 
and radius of the cylinder, respectively, r^s is the sth 
root of the Besscl hmction of order to, the normaliza- 
tion constant M depends on {m, s}, and n is an integer. 
In this basis set, the light-matter interaction matrix ele- 
ments read 



-iujt 



{ck'm'\h} ^ |wkTO) = 6cfe;,m',i;fc,m e 
{vk'm'l h\-^ |ckm) = 4\v™,.fe>^' e^"* x 

5-l,m'-rn5-q^,k'^-k^ , (H) 



where ^b'k'^m,bkrm = -;^Ao(p6'6 • e^)Af'Af drr x 

Jl{qr r)Ji+,n{K f) JmiK r). 

In order to derive the TL-generalized SBE, we special- 
ize Eqs. ([MT]) to cylindrical states, for which the compos- 
ite quantum index is a = {fc^, fc^, to}. The plan is to write 
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the equations of motion only for those components of p 
which are effectively coupled among themselves. In this 
way, we separate the evolution on the whole Hilbert space 
into that of dynamically uncoupled subspaces. Thus, we 
proceed by first writing down the equation of motion of a 
population or an intraband coherence, say, of the valence 
band. It can be seen that if the semiconductor is initially 



in its non-interacting ground state (it is universally ac- 
cepted in the SBE literature to assume the absence of 
Coulomb correlations in the unexcited material) we only 
need to consider cases with {fc^ = kz,m' = m}, since 
all other coherences remain zero at all times Thus we 
write 



7, Pvk^k' m. k.kr-m — ^vk^kj.7n, k-k^ m Pvk^k' m. k. kr-m ~^ 
at ' -r- 

iu>t \ ^ C* n e^i^i 

^ / J 'ick'Jm+l. vk^m P'uk-k'^m, ck^+q^k'^m+l e 



^ ^ ^ck'J m+l^ vkr-m Pck^+q^k'^m+l, vk^k^-m • (^^) 
fc" 



We see that Pvk,k'^m,k,k^m gets coupled to an interband coherence Pvk,k'^m,ck,+q,k';m+i and Pck-^+q,k';m+i, vk:,krm. with 
all values fc" of the radial quantum number, but with just {fc^,m} and {fc^ + (j^,to + for the other two quantum 
numbers. The Heisenberg equation for these interband coherences is 

^'^ ^^Pvkzk'^m, ckz+qzkrm+l ^ckz+qzk-r^n+l.vkzk'^ra Pvkzk'^ni, ckz+qzkrni+l ~t~ 

^ ^ ^ ^ckj.ra+1^ vk'^m Pvkzk^^m, kzk'^m ^ ^ck'^jn-i-l^ vk'^m Pckz+qzk'^ m+l, kz-hqzkrm-\-l j (-^^) 

k'^ k'^ 

with pva',ca = Pea va' ■ Inspection of this equation reveals that the interband coherence is coupled back to the initial 
valence-band coherence [Eq. (|12p ] and additionally to a conduction-band coherence or population, whose equation of 
motion is 

'^'^ ^^Pckz+qzk'^ra+l.kz+qzkrm+l ^ckz+qzkrm-{-l,kz-\-qzk'^7n+l Pckz+qzk'^7n-\-l,kz-\-qzkrm-\-l 

6 ^ ^^ S,ckz+qzk'^m+l,vkzk'Jm Pckz+qzk'^ni+l,vkzk'Jm ~ 6 ^ ^j ^ck^+qz.kr-rn+l.vk^k'Sm Pvkzk'^m,ckz+qzkrm+l ■ (14) 



This equation couples the population or intraband co- 
herence pckz+qzk'^m+i,kz+qzk^m+i to iutcrbaud cohereuces 
which evolve according to Eq. ([T^. It is clear that the 
system of equations is closed in the subspaces of fixed 
{kz,m} and {fc^ + qz , m + 1} , and the complexity of the 
problem has been drastically reduced, compared to the 
system of equations found in the case of Bloch states for 
a cubic bulk material. We may say that our procedure is 
equivalent to a block diagonalization. At this stage, the 
problem is highly tractable by computational techniques, 
since the only unconstrained variable is fc^.. 



C. Low-excitation regime 

Equations ([T^ - (|14p in all their generality are not 
amenable to analytical treatment. However, in the case 
of low photo-excitation, an analytical perturbative ap- 
proach is possible and gives us the basic physical insight 
that we are looking for. We now pursue this approach, 
but work initially with the system of equations (H])-© 
instead of ((T^ - (|14p . since the former are more general 



and also notationally simpler than the latter. We solve 
the system to lowest order in the vector potential, that 
is, we first solve Eq. (I7|) assuming that the zeroth-order 
intraband elements are Pv,a'a = Sa',a and Pc.a'a ~ 0, 
and then solve Eqs. ([5]) and (|6]) using the first-order so- 
lution of ([7|). The equation of motion for the interband 
polarization Pva',ca becomes 



ih- 



dt 



p'''\ 

I va' .CO. 



{ca\ hi \va') . (15) 



For a monochromatic electromagnetic field turned on at 
t = the solution reads 



pi'i'.caW = rca,™'(0(ca|4^^|W), (16) 



with 



^col.vol' (^) 



,)-hu] t/h 
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Inserting p^"^' ca ^^'^ equations for the intra-band coherence, Eqs. ([5]) and ([6]), we obtain 



')*/'^^(ca| hy' \vai){vai\hY> \ca') 

ai 

dt ""^ ^ ^ \Xca' ,vai{^ ) ~ Yca.vaiii )] 7 



dt "° ) / [ycQi,ua'(t ) — ^cQi.tral^ )] • 



For example, the conduction-band populations reduce to 



rc.a 



fc/ / ' 'J' {l~cos £,„-£.^, -MV^} 



(17) 

(18) 
(19) 

(20) 



where one sees that their time evolution is slow, with frequencies close to the detuning. As for the products of matrix 
elements, in our problem, with a = {fcz, fc^, m}, a simple calculation using Eqs. (jlip yields 



(ca| hY'^ \vOL-^{va-Y\'h\ ^ |ca') = Ccfc,m,Dfci,miCfc'm,i;fci,m/m',m4,,fc;,'5mi,m-i4i,,fc,-g, 



(wa| h\ ' |cai)(cai| = Cfei,mi,«fc,m?cfci,mi,i,fc; 



(21) 



r 



in agreement with the decoupling of Eqs. (|T2 l) -([T4 l) . show- 
ing that the second-order process involves an intermedi- 
ate state ai which can only differ from the initial state 
by ±Z and zbg^ in the quantum numbers m and A:^, re- 
spectively. Note that the interband coherence is of order 
0[A(r, i)^], while the populations and intraband coher- 
ences are of order C'[A(r, i)^], as indicated with super- 
scripts. Finally, the time behavior of each component is 
clearly discernible: while the interband coherence oscil- 
lates at the frequency of the TL field, the populations and 
intraband coherences do it typically at terahertz frequen- 
cies associated with interband Rabi flops and intraband 
energy differences. 



IV. ELECTRON QUANTUM KINETICS 



The solutions to Eqs. ([T2|) - ([T4| are the building blocks 
for constructing mean values of observables of interest. 
In the standard theory of optical transitions, where the 
light is assumed to be a plane wave and the dipole ap- 
proximation is made, once the time- and momentum- 
dependent density matrix is obtained, one calculates the 
macroscopic optical polarization and from it, for exam- 
ple, the optical susceptibility.— Under those assump- 
tions, the macroscopic polarization is just a spatially uni- 
form, time-dependent, function. 

By contrast, if the inhomogeneities of the field are 
taken into account (e.g. finite beam waist and oscillatory 
dependence in the propagation direction), the electronic 



variables acquire an intricate space dependence. The ex- 
citation of solids by TL beams also produces a space- 
dependent carrier kinetics which requires local variables 
for its description. In order to visualize the pattern of 
motion of the photo-excited electrons, we calculate in 
this Section the spatially inhomogeneous electric current 
density. Another useful variable, the transferred angular 
momentum, is instead a global magnitude that charac- 
terizes the TL-material interaction. Here we calculate 
their dynamics up to second order in the field amplitude. 

In the calculations that follow, we will study separately 
the contributions to the angular momentum and the elec- 
tric current made by the interband coherences, on the 
one hand, and by the populations and intraband coher- 
ences, on the other. This separation is conceptually use- 
ful because the interband coherence contributions are fast 
(femtosecond) oscillations around a null value of the cur- 
rent or angular momentum, analogously to what happens 
with the interband polarization, while the population or 
intraband-coherence contributions come from slower (pi- 
cosecond) processes in which a net transfer of momentum 
from light to matter occurs. The latter are related to the 
photon-drag effect fi^ii^ which is now generalized to in- 
corporate a rotational drag in the plane perpendicular to 
the propagation direction, due to the "slow" transfer of 
angular momentum. Furthermore, as we will see below, 
the lowest-order contributions for interband and intra- 
band processes are of first and second order in the light 
field, respectively. 

An investigation of the transfer of orbital angular mo- 
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mentum and the generation of paramagnetic currents in 
semiconductors and insulators was presented by us in 
a previous publication.— In that study we employed a 
simple wave-function approach, which was limited to de- 
scribe the single-particle dynamics. Our current the- 
oretical analysis uses the formalism of the Heisenberg 
equations of motion for populations and coherences, 
which fully accounts for the Pauli exclusion in the multi- 
electron excitation process, and has the advantage that it 
can be extended to include electron-electron interaction. 
In what follows, we analyze with this tool the transfer of 
orbital angular mometum and the generation of electric 
currents, including the diamagnetic term — missing in our 
previous study. 



A. Transfer of angular momentum 

Since the TL beam carries angular momentum around 
the z-axis. we focus on the corresponding quantity for 
electrons, 

= {b'a'\i,\ba)al,^,{t)ato,it), (22) 

b' haoi' 

where Iz = —ihd^. We split the matrix-element inte- 
gral {b'a'l Iz \ba) into an integral on the unit cell and a 
sum over lattice sites. Naturally, care must be taken 
when operating with —ihd^ on the envelope, 4>km(r) = 
MJmikr r)e*™'^e*'^^^, and on the periodic, u&(r), parts of 
the wave function ■0bkm(i")- We obtain 



(&'k'm'|;^|6km) 



(1 — (5b&') •^b'k'm',bkm , 



where 



(23) 



and this integral is over the whole crystal. We note that 
the quantity lz,b'b does not depend on the 0AM of light, 
and so we disregard it from now on. Next, we split the 
angular momentum into interband (coherence) and in- 
traband (population and coherence) contributions, and 
use the perturbation-theory solutions of Eqs. (fT6l) -(fT9 |) . 

Interband- coherence contribution: After taking the 
mean value of the angular momentum operator [written 
in second-quantization notation in Eq. ([22])] on the initial 
state, we identify the component of the electronic angular 
momentum driven by the interband coherence as 



k'rn'km 



23? 



■C-uk'm'.ckm P^k'm',ck;m(^) -(24) 



We note that £uk'm',ckm contains the factor Ski^k^, while 
Pvk'7n',cUm(t) contains the factor Sk,,k'^+q,. This mis- 
match, which renders a vanishing L 



(coh) 



comes from 



dropping the dipole approximation in the calculation of 
the interband coherence, and, at the same time, consid- 
ering that the solid is infinite in extent [in Eq. ([23)) ]. To 
be consistent, one needs to consider the system as a thin 
slice of semiconductor perpendicular to the z-axis, having 
a width much smaller than the wavelength of the light. 
As a consequence, if the slice is located at zq, Sk',^k, 
replaced by exp[i(kz — k'^)zo] and we obtain 



IS 



^b'Wni' 



m' .bkm 



i{k^-ki)zo y 
k^m e A 



' rn' ,kj.ri\ 



i TT fk 

[5-l,m-m'P-,b'b — Si, 



'P+,b'b\ , 



with fk'^rn',krm = M' M J dr r"^ Jm' (K r) JmiK r) , and 
P±,b'b = X ■ pb'b ± • Pb'b- Then, we may succinctly 
state that /Cb'k'm'.bkm S±i,m-m'- luspcctiou of the co- 
herence, Eq. ()16p . and matrix elements, Eq. shows 



Thus, wc conclude that at 



this level there is a transfer of angular momentum back 
and forth between the light beam and the electrons if and 
only if |Z| = 1. We emphasize that, on time-average, there 
is no net transfer of angular momentum to the material 
system, unless the temporal shape of the electromagnetic 
pulse is asymmetric 1^ 

Population and intraband- coherence contribution: The 
component of angular momentum driven by the popula- 
tions and intraband coherences reads 



L'f°P\t) 



km 



hm 



(2) 
-ukn 



(2) 
Pckn 



.it) 



In order to correctly interpret this expression, we re- 
call that the TL photo-excitation process connects 
the valcncc-band subspacc of fixed {kz,m} with the 
conduction-band subspacc of fixed {fc^ + Qz, m + 1}; thus, 
an imbalance population in the conduction band is pro- 
duced. This asymmetry between populations in both 
bands with respect to the quantum number m brings 
about a net angular momentum acquired by the elec- 
trons, which we refer to as rotational photon-drag. In 



contrast to L^z°'^\t), Lf"'''{t) has no restrictions on the 
values of I that cause a transfer of angular momentum, 
and its time average yields a non-zero value. 



- (pop) 



B. Induced currents 



Next wc will obtain the photo-induced currents pro- 
duced by the irradiation with TL. The general ex- 
pressions of the electric current in second quanti- 
zation notation are as follows. The more stan- 
dard, paramagnetic, current density is given by 
j(p)(x,t) = -i(7ft/(2me) W^,(V- V')^/'t(a;',t)7/;(a;,t) 



and the diamagnetic term is given by 'y >{x,t) = 
—q/meA{r,t)tp''{x,t)i(j{x,t)^. We apply these expres- 
sions to our problem, and after some algebraic manipu- 
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lation, wc obtain, for the paramagnetic term: 



Interband- coherence contribution: The interband- 
coherence contribution to the current density is given by 



2 nip 



6km 



rr) X 



V-bkmW VV-b'k'm'W] ai'k'm'(Oabkm(i), (25) 
and for the diamagnetic term: 

977 _ ' ^ 



b'k'm' 
fekm 



«h'k'm'(0 a(,km(<) . 



(26) 



In what follows both contributions will be studied in de- 
tail. The simplification of the expressions will proceed 
in a similar manner to the calculation of the transferred 
angular momentum. However, at a certain point we will 
make use of a space average (A = (l/a"^) J^^^ d'^r A) in 
order to eliminate irrelevant microscopic (intra-cell) de- 
tails. 



1. Paramagnetic- current density 



Now we work out the expression psp by replacing the 
wave functions from Eq. (fTO|). We apply the gradient op- 
erator on the envelope and periodic parts of the Bloch 
wave function, perform space average, take mean value 
over the initial state, and split the result into interband- 
coherence and population and intraband-coherence con- 
tributions. 



Y h 



X 



FIG. 1: (color online) First-order paramagnetic current for 
light with / = 1 and polarization a — . The center of the plot 
coincides with the beam axis. The factor Jm-i{k'j.r)Jm.{krr) 
is drawn in dashed (red) line. 



k'm' 
km 

i{k^-k'^)z i{m-m')4> {i) 
^ ^ f^uk'm'.ckmV'''' 



With the help of Eqs. (fTTj) and (|16p we simplify this ex- 
pression to 



-{coh) 2(7 

J (r^O — ^ 



(1) 



Jm~l{k^r)Jm{krr) Pyk^-q^k'^m~l,ck^k^mi^) 



. . / ^ ^ \ V . 

' q.z=0 ^ 



V \ i i / . 



' q.z=7: ' 



\ 



\ 



(\.z=nl2 



^q.z=3;r/2' 



FIG. 2: (color online) First-order paramagnetic current for 
light with / = 1. The center of each plot coincides with the 
beam axis. Each panel pictures a different slice along the 
2-axis, showing the wave nature and space variation of the 
pattern. 

The main feature of this expression is that it contains 
a sum over products of space- and time-dependent func- 
tions. The time dependence presents two distinct scales, 
as mentioned before: a rapid oscillation related to the 
frequency uj of the light beam, and a slower one that is 
related to the detuning. 

To describe the electric-current patterns in the plane 
perpendicular to the z-axis, we disregard the slow time 
evolution and focus on the space- and time-dependent 
quantity 5R[p,ee^(«^"+''^-"*) ^^-((fc^r) J„(fc,r)], see Fig. 
[T] We can see that the current density and the angular 
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momentum at this level are consistent with each other. 
To see this, let us take as an example a TL field having 
circular polarization cr-|_ = x + iy. From the usual selec- 
tion rule for the absorption of a photon, a non-vanishing 
light- matter matrix element requires = po{x — iy). 
Since p„c = P™ ^ Poi^ + iy), we see from Eq. ([24|) that 
a non-vanishing angular momentum will appear only if 
Z = — 1; in this case the current pattern reflects this fact, 
presenting a "circulation" around the beam axis. On the 
contrary, if the beam is tuned to ^ — 4-1 and a+, the cur- 
rent forms a pattern that docs not "flow" around the axis. 
As we move along the z-axis for fixed time, we observe 
the wave nature of the x-y plane current (Fig. [2]). For 
fixed z, the wave nature is revealed as time evolves. For 
different values of I, the electric current develops more 
complex patterns, which mimic the complex structure of 
the electric field of a TL beam. In the case of |^| > 1, 
there appear more than one off-centered vortices, as il- 
lustrated in Fig. m 



X \ Y Ji / / ^ 



paramagnetic current is given by 



/ ■ 

/ - 
\ 

\ - 



I i 1 1 ' \ \ l \ — f 




/ f f t \ \ \ 



FIG. 3: (color online) First-order paramagnetic current for 
light with 1 = 2. The center of the plot coincides with the 
beam axis. The space pattern is more complex than in the 
I = 1 case and exhibits circulation around off-axis centers. 



Prom the above we can extract some general features 

of i^'^°'^\r,t). Driven by the interband coherence, it os- 
cillates in time, with zero mean, at the frequency of the 
EM field, and presents complex spatial patterns that are 
related to the peculiar electric field of the TL beam. For 
special cases, the spatial pattern displays a circulation (of 
microscopic origin and not to be confused with a macro- 
scopic excursion of the electrons around the beam axis) 
related to the non-vanishing coherence contribution of 
angular momentum calculated in Sec. IIV Al 

Population and intraband- coherence contribution: The 
population and intraband-coherence contribution to the 



-{pop} I ,^ 



q h 
2 mp 



J2 {[$k'™'(r)V$k™(r)]- 



k'm' 
km 

(2) 



[$k™(r)V<I>J:,,„,(r)]} plZ,n'M,nii) + {c^v}, 

where PcWm'.kmW is given by Eqs. (HD) and ([21]), and 
{c — > w} stands for a similar term replacing c by w. The 
intraband current in the direction of 6 is 



Af'Afm-Jni(kLr)Jm(krr) 



rakz 

(2) 

'^ck'kzlTl-.krk- 771 



(t) + {c^v}. 



(27) 



Given that the electrons excited by the TL beam occupy 
initially a portion of the valence band that is symmetric 
with respect to the F point of the Brillouin zone, we have 

disregarded the contribution to (r, t) coming from 

the holes left behind, and kept only the current produced 
by the imbalance of electrons in the conduction band. 
The parameter I does not appear explicitly, but it enters 

(2) 

Pck' k-m krk,mi^)^ sincc an electron leaving a state with 
{v, m} goes to a state with {c, m -\- 1}. 



2. Diamagnetic current density 

Starting from Eq. (|26p . we perform a space average 
and obtain for the diamagnetic term 



j('^)(r,t) 



A(r,t) ^f'^fJra{K^)Jm{krr) 



k'^k^ 
k~m 



PcKk,m,k^k,mii) + {C ^ ^ 

We point out the following features of the diamagnetic 
current density: i) it arises from populations and intra- 
band coherences in each band; ii) its vectorial character 
is given by the polarization of the light; Hi) it is of third 
order in the vector potential amplitude; iv) it does not 
arise from an imbalance of the conduction-band popula- 
tion [it is not proportional to m like j^^°^\v,t)^ see Eq. 
(|27l) ]: v) its time evolution is given by the slow evolution 

(2) 

of pc (t) and the fast oscillation of the field A{r,t). 



V. CONCLUSIONS 

We developed a theory of the band-to-band transitions 
induced by twisted light (light carrying orbital angular 
momentum) in bulk semiconductors. Wc posed the prob- 
lem of the light-matter interaction in terms of Hcisenberg 
equations of motion for the populations and quantum 
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coherences of a two-band semiconductor model, as cus- 
tomarily done. We found that the resulting system of 
equations is greatly simplified when the envelope electron 
wave functions are represented in cylindrical coordinates, 
instead of using the usual (cartesian) Bloch state repre- 
sentation. This simplification is due to the decoupling 
of the system of equations in subsystems determined by 
the orbital angular momentum of electrons. Though non- 
standard for the bulk case, our choice of basis states is, on 
the one hand, perfectly admissible and, on the other, it 
proves to be the best choice from a mathematical point of 
view. It is admissible because the material properties are 
unaffected by surface effects in the limit of a bulk/large 
system. It is the right choice since it provides the highest 
symmetry compatibility between the twisted-light vector 
potential and the electron states. 

Despite the achieved simplification, the evolution of 
the different relevant physical quantities under the exci- 
tation by a time-dependent pulse must be computed via 
numerical analysis of the equations of motion. This task 
is left for future work; instead, we showed here analyti- 
cal results in the low-excitation or pcrturbative regime. 
With the solutions for the populations and quantum co- 
herences, we confirmed on more solid grounds our pre- 
vious findings, i.e. that the optical excitation will gener- 
ate electric currents, and that there will be a transfer of 
orbital angular momentum from the light beam to the 
electrons. Our analysis of the electric current and elec- 
tron's orbital angular momentum showed that two qual- 
itatively different contributions enter both observables; 
they may be termed microscopic and macroscopic con- 
tributions. The microscopic contribution relates to the 
interband coherence and mimics the behavior of the elec- 
tric field; from this and other reasons, it parallels the in- 
duced polarization of a semiconductor in the presence of 
plane- waves, as traditionally studied using the vertical- 
transition assumption. On the other hand, the macro- 
scopic contribution signals a net transfer of 0AM from 
the field to the electrons, and it parallels the photon- 
drag effect. We showed that the electric current exhibits 
a high degree of spatial complexity, due to the inhomo- 
geneous nature of the twisted-light beam. Additionally, 
we have calculated and briefiy analyzed the diamagnetic 
current, not addressed in our previous study. 



known solution of the Schrodinger equation without po- 
tential 



2m 



V2$(r) = E^{r) , 



where the Laplacian in cylindrical coordinates is 



ld_ 

r dr 



,dl 
dr 



1 d^f d^f 

J.2 QQ2 Q^2 



Consider a solution of the separable form <I'(r) = 
R{r)Q{0)Z{z); replacing into the Schrodinger equation, 
we get 



2m 



1 1 d 
Rr dr 



OR 
dr 



1 1 a^e 1 d'^z 



e r2 9612 z 9z2 



which yields 



" 1 1 


d 




1 d^Z' 


1 1 d^Q 




dr 




^ Z dz\ 


+ e ,,2 QQ2 



= E 



-a 



where = a^. Then 



1 d ( dR\ r2 d'^Z 



— r — r — 

R dr \ dr 



Z d£^ 



2 2 

+ Q! r 



1 929 
e 9612 



which splits to 



1 A 
'R'd^^ 



dR 
dr 



Z dz^ 



22 2 

a r = m 



1 929 
6 9^ 



The second equation has solution 9 = Aie""^ -1-^26 
The remaining equation is 



-iraO 



l_ld_ 

R r dr 



dR 
dr 



1 d^z 

Z'd^ 



= 0. 



This is again separable 



1 1 9 

Rr dr 



dR 



2 ^2 

T ~ 



1 d^Z 
Z'd^ 
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with solution Z = Be , and equation 



I d f dR 



— r — r 

R dr \ dr 



+ [a^ - - m^ 



0, 



and developing the derivatives 
^d?R dR 



dr^ 



dr 



+ [{a^ ~ X^y - m"^] R = 0, 



Appendix A: Particle in a hollow cylinder 



which is almost the Bessel differential equation. Defining 
^2 = (a2 _ A2)r2 the correct differential equation follows 



We present the complete derivation of the electronic 
states in cylindrical coordinates, starting from the well- 



^d^R , dR 



da; 2 



x-P + (x^ - m^) i? = 0. 
dx ^ ' 
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1. Boundary conditions 

For Z{z) we use boundary conditions Z{Q) — Z{H), 
and obtain 



where U (r) is the lattice potential Since we have already 
solved the problem of the free particle without u{r), we 

know that -^V2$„ms(r) E^Xs^nmsir)- Dividing 

by ^r, 



(r) and grouping terms we get 



Z(z) ^ ^e'^'- 



X = 



2Tm 

IT 



For 6(6'), the periodicity implies 

6(61) = iVee'""* 

m = ...,-2,-1,0,1,2,... . 

For the radial solution, since it has to be finite at the 
origin, the solutions are the Bessel functions of the first 
kind Jm{x)- If we demand that Jm{x = ^ — }?V) = 0, 
then the sth zero rms of the Bessel function should be 

^2 ' - 



En 



2m 



27rn 

IT 



,(r) = MnsJ,„(V«2-A2 r)e™^e' 



with normalization Afms = 
for any admissible solution. 



and A < 



2. Cylinder with a Bravais lattice 

In the effective-mass approximation, the complete 
wave function is expressed as the product of an enve- 
lope $„ms(r) and a periodic u{r) function. Then, the 
Schrodinger equation reads 

-— it(r)V2*„™3(r) - — $„„^(r)V2u(r) - 
2m 2m 



— V$„™^ 
m 



(r)] • [Vu(r)] + C/(r)[$„„,(r)?/(r)] = 

E^{y)u{v) , 



-— V^u r + U{r)u(v) 

2m m 



,(r) • Vii(r) 



^nms(r) 



[E-E^lM^)- 



This is an equation for u{v) [since we already know the 
functional form of $„ms(r)]. We ask that u{t) — u{r + 
R), with R a lattice vector, so we can regard the region of 
integration as the unit cell. Then, if $nms varies slowly 
in the unit cell, we see (?i^/m)<J>„„is(r)~^ [V$„„is(r)] • 
[VM(r)] as a perturbation. This is simply the k • p ap- 
proximation in a different coordinate system. To lowest 
order we have 



2m 



V^u{r) + U{r)u{r) = AEu{r). 



and from here one obtains u(r) as usual. Therefore, we 
take the eigenfunctions as a product of envelope and pe- 
riodic functions, while the energy is that of the envelope 
corrected by the energy gap and the effective mass: 



E, 



brims 



2ml 



SbcEg 



Nr 



where H is the height and L the radius of the cylinder. 
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